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Since the detection of the asymptotic properties of the dipole gravity modes in the Sun, the quest to find the individual 

gravity modes has continued. A deeper analysis of the GOLF/SoHO data unveils the presence of a pattern of peaks 

that could be interpreted as individual dipole gravity modes. The computed collapsed spectrum -around these candidate 

modes- uncovers the presence of a quasi constant frequency splitting, in contrast with regions where no g modes are 

expected in which the collapsogram gives random results. Besides, the same technique applied to VIRGO/SoHO unveils 

some common signals between both power spectra. Thus, we can identify and characterize the modes, for example, with 

their central frequency and splittings. This would open the path towards new investigations to better constrain the solar 
core. 
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1 Introduction 

Our knowledge of the Sun's interior has been considerably 
improved by the observations of the solar acoustic oscilla- 
tions (p modes). Indeed, accurate measurements of the low- 
degree p modes (Toutain et al. 1997; Thiery et al. 2000; 
Garcia et al. 2001; Salabert et al. 2004; Broomhall et al. 
2009) have allowed us to infer the structure (e.g. Turck- 
Chieze et al. 2001; Basu et al. 2009) and dynamics (Chaplin 
et al. 2001; Garcia et al. 2004; Mathis & Zahn 2004, 2005; 
Turck-Chieze et al. 2010) of the layers below the convec- 
tive zone, from the radiative zone down to the solar core. 
However, the Sun's deepest regions are far from being fully 
known, and the inner radiative zone below 0.25 i? , where 
about 50% of the total mass is concentrated, will be only in- 
ferred by the detection of individual gravity (g) modes. For 
instance, few g modes would considerably improve the in- 
versions of the rotational profile down to 0. 1 Rq (Mathur et 
al. 2008, 2009), while we also need to improve the detection 
of high-order low-degree p modes (Garcia et al. 2008c). The 
increasing background convective level towards the lower 
frequencies (e.g. Lefebvre et al. 2008), combined with very 
small amplitudes of those modes (several mm/s in the case 
of g modes, Belkacem et al. 2009) are the limiting fac- 
tors for their detections. In the case of the low-degree, low- 
frequency p modes, accurate measurements are hardly ob- 
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tained below 1 mHz (Bertello et al. 2000; Broomhall et al. 
2008, 2010; Salabert et al. 2009), while no general con- 
sensus has been obtained for the detection of individual g 
modes (Appourchaux et al. 2010). It is necessary then to im- 
prove the current analysis methods and to develop new ones 
by trying to increase the signal-to-noise ratio (SNR) at low 
frequency. That can be achieved, for example, by reducing 
the convective noise (e.g. Garcia et al. 1999), by combining 
several instruments (e.g. Broomhall et al. 2007, 2010; Sal- 
abert & Garcia 2008), or by using new instruments, such as 
PICARD (Thuillier et al. 2006), GOLF-NG (Turck- Chieze 
et al. 2006) and the very promising HMI and AIA aboard 
SDO (Howe 2010). Nevertheless, signals in the g-mode re- 
gion have been detected in the GOLF and VIRGO data on- 
board SoHO, which are very unlikely to be due to instru- 
ment or solar noise (Jimenez & Garcia 2009). Moreover, 
several patterns with high confidence level have been ob- 
served above the solar convective noise (Turck- Chieze et 
al. 2004; Garcia et al. 2008a). The observation of the pe- 
riodic signature of the g modes can be obtained by using 
methodologies applied in asteroseismology, where the large 
separations of solar-like stars with very small SNR can be 
measured even without being able to see clear features in 
the power spectrum (e.g. Mosser et al. 2009; Mathur et al. 
2010). This technique has provided the first detection of the 
asymptotic periodicity of dipole g modes by computing the 
power spectrum of the power spectrum (PSPS) (Garcia et 
al. 2007). From this analysis, it was possible to infer an av- 
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erage speed in the solar core (3 to 5 times faster than the rest 
of the radiative zone), and an internal structure compatible 
with solar models based on the "old" abundances (Garcia et 
al. 2008b). 

2 Observations 

We have used 4472 days of Global Oscillations at Low Fre- 
quency (GOLF; Gabriel et al. 1995) time series calibrated 
into velocity (Palle et al. 1999; Ulrich et al. 2000; Garcia et 
al. 2005), starting April 11, 1996. We have worked with a 
single - full resolution - power spectrum even knowing that 
GOLF has been observing in two different configurations 
with a different sensitivity to the visible solar disk (Garcia 
et al. 1998; Henney et al. 1999). To reduce the problem of 
the discretization in frequency, we have computed a 5 times 
zero-padded power spectrum (see the discussion concern- 
ing this problem in Gabriel et al. 2002). We have also used 
time series from the SPM/VIRGO package (Frohlich et al. 
1995). 
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Fig. 1 5 times zero-padded PSD of 4472 days of GOLF 
velocity data in which a smooth by a boxcar of 41 nHz has 
been done. The vertical red dotted lines are the 1=1 cen- 
tral (m=0) g-mode frequencies computed using the Saclay 
seismic model (Mathuret al. 2007). Dashed and dot-dashed 
vertical red lines are the m=l components assuming a core 
rotating 4.5 times the rest of the radiative region. 



3 Data analysis 

In 2007, Garcia et al (2007) uncovered the existence of a 
pattern of peaks in the power spectrum of the GOLF in- 
strument in the region between 25 and 140 /iHz that were 
equidistant in period. They interpreted this periodicity as 
the asymptotic properties of the dipole gravity modes. How- 
ever, they were not able to distinguish the individual modes 
in the power spectrum. To face this problem and to try to 
identify the individual 1=1 g modes we follow in this paper 
the same methodology commonly used in asteroseismology. 
Indeed, the first parameters that are searched for in astero- 
seismology are the frequency separations as p modes are 
equidistant in frequency (e.g. Bedding et al. 2010; Stello et 
al. 2010), while when the SNR is high enough, we can prop- 
erly determine the characteristics of the individual p modes 
(e.g. Arentoft et al. 2008; Appourchaux et al. 2008; Garcia 
et al. 2009; Chaplin et al. 2010; Deheuvels et al. 2010). On 
the other hand, when the SNR of the modes is very small we 
are still able to determine the large separation even with- 
out seeing the individual p modes in the power spectrum 
(e.g. Mosser et al. 2009; Mathur et al. 2010). In such cases, 
the only thing that is generally done is to perform a heavy 
smooth of the power spectrum to unveil the p-mode hump 
in the region in which the frequency spacing is found. 

Figure [T] shows the 5 times zero-padded Power Spec- 
tral Density (PSD) of the GOLF time series. To increase 
the SNR of the possible signals we have applied a smooth 
of the spectrum by a boxcar function of 41 nHz. To guide 
the eyes, we have superimposed the predicted frequencies 
of the 1=1, m=0 components of the dipole gravity modes 
(vertical red dotted lines) from the Saclay seismic model 
(Couvidat et al. 2003; Mathur et al. 2007). In this frequency 
range, the differences between this solar model, the model S 
(Christensen-Dalsgaard et al. 1996) and the Ml model from 



Nice (Provost et al. 2000) are less than 0.05 ^Hz -which is 
inside the width of the vertical lines- and we have also veri- 
fied that the sensitivity of the pulsation codes to their param- 
eters is negligible for the 1=1 modes in this frequency range 
(Moya, Mathur, Garcia 2010). The dashed and dot-dashed 
vertical red lines are the m=±l assuming a core rotating 4.5 
times faster than the rest of the radiative region. This value 
is compatible with the range of the core rotation inferred by 
Garcia et al. (2007) and it matches the highest peaks of the 
PSD. Indeed, If we filter out all the peaks in such a way that 
we keep only the highest 0.3% of the peaks in the PSD and 
then we also remove the remaining spikes that are further 
than ±2.5 /iHz of the predicted 1=1 frequencies, we obtain 
the PSD showed in Fig. [2] 
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Fig. 2 Same than Fig. 1 but after filtering out all the peaks 
in the PSD in such a way that we only keep the highest 0.3% 
of the peaks in the PSD and also keeping the remaining 
peaks around ± 2.5 /iHz of the 1=1 theoretical frequencies. 
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Now, we can clearly identify in the PSD the pattern of 
peaks -equidistant in period- that was responsible of the de- 
tection of the APi spacing by Garcia et al. (2007) with more 
than 99.99% confidence level. Moreover, in 6 out of the 7 
modes analyzed in this frequency region, the peaks seems to 
be split in frequency which is what we would have expected 
for gravity modes traveling in a region where the internal 
rotation is high. To perform a more quantitative analysis of 
the frequency splittings we compute the collapsograms. 

4 Collapsograms 

The so called m-collapsograms were originally developed 
for imaged instruments. This technique consists of shifting 
each m-component of the PSD by a quantity proportional 
to the splitting and co-adding all of them. The resultant av- 
erage spectrum is fitted assuming a Lorentzian profile. The 
fitted peak with the smaller linewidth or maximum likeli- 
hood is the best result, which corresponds to a given dis- 
placement of the m-components: the searched splitting. The 
central frequency of the mode (m=0) is given by the fre- 
quency of the fitted Lorentzian (see Salabert et al. 2009 for 
a detailed explanation of the method). 

In the case of GOLF, we only have one PSD in which 
only the even (1 + m) components of the modes are visible. 
Thus, to mimic the methodology, we duplicate the GOLF 
PSD and we simply shift one of them in respect to the other 
one (like an auto correlation function). For each displace- 
ment we add both spectra and we look for the highest peak. 
The shift giving the highest peak in the collapsed spectrum 
provides the splitting and the frequency of the highest peak 
is the frequency of the mode. An example of the resultant 
collapsed spectrum is shown in Fig.[5]for the mode 1=1, n=- 
10. 



cases the maximum is obtained at the same frequency but 
with slightly different frequency shift (splitting) of less than 
10% (2 frequency bins). 
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Fig. 3 Collapsed spectrum around the mode 1=1, n=-10 
showing the displacement giving the highest peak (895 
nHz). The triple dot-dashed line marks the position of the 
predicted frequency. 

Figure|4]shows the collapsed spectrum for the mode 1= 1 , 
n=-8 for GOLF (top) and VIRGO (Bottom) data. In both 
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Fig. 4 Collapsed spectrum around the mode 1=1, n=-8 for 
GOLF (Top) and VIRGO/SPM (Bottom). 



Figure[5]shows the resultant frequency shifts (splittings) 
of the collapsograms of 6 out of 7 dipole modes in the region 
60 to 140 /iHz (black asterisks). A stable splitting around 
850-900 nHz has been found for all the modes. To do the 
collapsograms, we have scanned ± 1/iHz around the theoret- 
ical frequency, and rotation rates in a range 0.5 to 5 times the 
rotation of the rest of the radiative zone. Two modes were 
more difficult to characterize: the 1= 1, n= -7 ( 84/iHz) and 
the n= -5 ( 109 //Hz) because the m=-l and m=+l compo- 
nents are not clearly visible simultaneously. For the second 
mode (n=-5), we managed to obtain a splitting in the same 
range than the other modes by taking the second highest 
peak in the collapsogram. The blue diamonds in Fig. [5] rep- 
resents the collapsograms computed in regions in between 
the 1=1 modes where only noise or high-degree g modes 
would be present. A high dispersion of the results is found 
in the frequency shifts, which is what we expected for re- 
gions dominated by noise. 

5 Conclusions 

For the first time we have been able to identify and char- 
acterize the individual peaks responsible for the detection 
of the asymptotic APi spacing measured by Garcia et al. 
(2007) with more than 99.99% confidence level. Moreover, 
the peaks showed a quasi constant frequency splitting in a 
range 850 to 900 nHz. This work is still in progress, but if 
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Fig. 5 Splittings obtained for the 6 out of 7 1= 1 g modes in 
the regin 60-140 /iHz (Black asterisks). The blue diamonds 
ae the shifts obtained in regions where we do not expect 
gravity modes. 



these results are confirmed, it would allow us to better study 
the core of the Sun 

Acknowledgements. GOLF and VIRGO instruments on board the 
SoHO spacecraft are a cooperative effort of many individuals, to 
whom we are indebted. SoHO is a space mission of international 
cooperation between ESA and NASA. D.S. acknowledges the sup- 
port of the grant PNAyA 2007-62650 from the Spanish National 
Research Plan, A.M. the support from a Juan de la Cierva contract 
of the Spanish Ministry of Science and Innovation and J.B. the 
support through the ANR SIROCO. R.A.G. and S.T thanks the 
support form CNES. 

References 

Appourchaux, T, Michel, E., Auvergne, M., et al.: 2008, A&A 

488, 705 
Appourchaux, T, Belkacem, K., Broomhall, A.M., et al.: 2010, 

ARA&A 18, 197 
Arentoft, T, Kjeldsen, H., Bedding, T, et al.: 2008, ApJ 687, 1 180 
Basu, S., Chaplin, W. J., Elsworth, Y., New, R., Serenelli, A. M.: 

2009, ApJ 699, 1403 
Bedding, T.R., Huber, D., Stello, D., et al.: 2010, ApJ 713, L176 
Belkacem, K., Samadi, R., Goupil, M. J., et al.: 2009, A&A 494, 

191 
Bertello, L., Varadi, E, Ulrich, R. K., et al.: 2000, ApJ 537, L143 
Broomhall, A. M., Chaplin, W. J., Elsworth, Y., Appourchaux, T: 

2007, MNRAS 379, 2 

Broomhall, A. M., Chaplin, W. J., Elsworth, Y, Fletcher, S. T: 

2008, Astronomische Nachrichten 329, 461 

Broomhall, A.-M., Chaplin, W. J., Davies, G R., et al.: 2009, MN- 
RAS 396, L100 

Broomhall, A.M., Chaplin, W.J., Elsworth, Y, et al.: 2010, MN- 
RAS, in press 

Chaplin, W. J., Elsworth, Y, Isaak, G R., et al.: 2001, MN- 
RAS 327,1127 

Chaplin, W.J., Appourchaux, T, Elsworth, Y, et al.: 2010, 
ApJ 713, L169 

Christensen-Dalsgaard, J., Dappen, W., Ajukov, S. V, et al.: 1996, 
Science 272, 1286 



Couvidat, S., Turck-Chieze, S., Kosovichev, A.G.: 2003, ApJ 599, 

1434 
Deheuvels, S., Bruntt, H., Michel, E., et al.: 2010, A&A 515, A87 
Frohlich, C, Romero, J., Roth, H., et al.: 1995, Sol. Phys. 162, 101 
Gabriel, A.H., Grec, G, Charra, J., et al.: 1995, Sol. Phys. 162, 61 
Gabriel, A.H., Baudin, F, Boumier, P., et al.: 2002, A&A 390, 

1119 
Garcia, R.A., Roca Cortes, T, Regulo, C: 1998, A&AS 128, 389 
Garcia, R. A., Jefferies, S. M., Toner, C. G, Palle, P. L.: 1999, 

A&A346, L61 
Garcia, R. A., Regulo, C, Turck-Chieze, S., et al.: 2001, Sol. 

Phys. 200, 361 
Garcia, R. A., Corbard, T, Chaplin, W. J., et al.: 2004, Sol. 

Phys. 220, 269 
Garcia, R.A. Turck-Chieze, S., Boumier, P., et al.: 2005, A&A 442, 

385 
Garcia, R. A., Turck-Chieze, S., Jimenez-Reyes, S. J., Ballot, J., 

Palle, P. L., Eff-Darwich, A., Mathur, S., Provost, J.: 2007, 

Science 316, 1591 
Garcia, R. A., Jimenez, A., Mathur, S., et al.: 2008a, Astronomis- 
che Nachrichten 329, 476 
Garcia, R. A., Mathur, S., Ballot, J.: 2008b, Sol. Phys. 251, 135 
Garcia, R. A., Mathur, S., Ballot, J., et al.: 2008c, Sol. Phys. 251, 

119 
Garcia, R.A., Regulo, C, Samadi, R., et al.: 2009, A&A 506, 41 
Henney, C.J., Ulrich, R.K., Bertello, L„ et al.: 1999, A&A 348, 

627 
Howe, R.: 2010, Proc. of GONG 2010/SOHO XXIV, In press 
Jimenez, A., Garcia, R. A.: 2009, ApJS 184, 288 
Lefebvre, S., Garcia, R. A., Jimenez-Reyes, S. J., Turck-Chieze, 

S., Mathur, S.: 2008, A&A 490, 1143 
Mathur, S., Turck-Chieze, S., Couvidat, S., Garcia, R. A.: 2007, 

ApJ 668, 594 
Mathur, S., Eff-Darwich, A., Garcia, R. A., Turck-Chieze, S.: 

2008, A&A 484, 517 
Mathur, S., Garcia, R. A., Eff-Darwich, A.: 2009, ASSP, 

<arXiv:09024142) 
Mathur, S., Garcia, R.A., Catala, C, et al.: 2010, A&A In press 
Mathis, S„ Zahn, J.-P: 2004, A&A 425, 229 
Mathis, S., Zahn, J.-P: 2005, A&A 440, 653 
Mosser, B., Michel, E., Appourchaux, T, et al.: 2009, A&A, 506, 

33 
Moya, A., Mathur, S., Garcia, R.A.: 2010, Sol. Phys. In press 
Palle, PL., Regulo, C, Roca Cortes, T, et al.,: 1999, A&A 341, 

625 
Provost, J., Berthomieu, G, Morel, P.: 2000, A&A 353, 775 
Salabert, D., Fossat, E., Gelly, B., et al.: 2004, A&A, 413, 1135 
Salabert, D. & Garcia, R. A.: 2008, ArXiv e-prints, 

larXiv:0810. 16961 
Salabert, D., Leibacher, J., Appourchaux, T, & Hill, F: 2009, 

ApJ 696, 653 
Stello, D., Basu, S., Bruntt, H, et al.: 2010, ApJ 713, L182 
Thiery, S., Boumier, P., Gabriel, A. H., et al.: 2000, A&A 355, 743 
Thuillier, G, Dewitte, S., Schmutz, W.: 2006, in 36th COSPAR 

Scientific Assembly, 170 
Toutain, T, Appourchaux, T, Baudin, F, et al.: 1997, Sol. 

Phys. 175,311 
Turck-Chieze, S., Couvidat, S., Kosovichec, A.G., et al.: 2001, 

ApJ 555, L69 
Turck-Chieze, S., Garcia, R. A., Couvidat, S., et al.: 2004, 

ApJ 604, 455 
Turck-Chieze, S., Carton, P.-H., Ballot, J., et al.: 2006, Advances 

in Space Research 38, 1812 



© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Wcinhcil 



www.an-journal.org 



Astron. Nachr. / AN (2006) 793 

Turck-Chieze, S., Palacios, A., Marques, J.P., Nghiem, P.A.P.: 

2010, ApJ 715, 1539 
Ulrich, R. K., Garcia, R. A., Robillot, J. M. et al.: 2000, A&A 364, 

799 



WWW.an-J0urnal.org © 2006 WILEY- VCH Verlag GmbH &Co. KGaA, Weinheim 



